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ABSTRACT 

Hysteresis  in  smart  materials  hinders  the  wider  applicability  of  such  materials  in  actuators.  In  this  paper,  a  systematic 
approach  for  coping  with  hysteresis  is  presented.  The  method  is  illustrated  through  the  example  of  controlling  a 
commercially  available  magnetostrictive  actuator. 

We  utilize  the  low-dimensional  model  for  the  magnetostrictive  actuator  that  was  developed  in  earlier  work.  For 
low  frequency  inputs,  the  model  approximates  to  a  rate-independent  hysteresis  operator,  with  current  as  its  input 
and  magnetization  as  its  output.  Magnetostrictive  strain  is  proportional  to  the  square  of  the  magnetization.  In  this 
paper,  we  use  a  classical  Preisach  operator  for  the  rate-independent  hysteresis  operator. 

In  this  paper,  we  present  the  results  of  experiments  conducted  on  a  commercial  magnetostrictive  actuator,  the 
purpose  of  which  was  the  control  of  the  displacement/strain  output.  A  constrained  least-squares  algorithm  is  em¬ 
ployed  to  identify  a  discrete  approximation  to  the  Preisach  measure.  We  then  discuss  a  nonlinear  inversion  algorithm 
for  the  resulting  Preisach  operator,  based  on  the  theory  of  strictly-increasing  operators.  This  algorithm  yields  a 
control  input  signal  to  produce  a  desired  magnetostrictive  response.  The  effectiveness  of  the  inversion  scheme  is 
demonstrated  via  an  open-loop  trajectory  tracking  experiment. 

Keywords:  Control,  hysteresis,  Preisach,  identification,  inversion,  magnetostriction,  smart  actuator 

1.  INTRODUCTION 

Hysteresis  in  smart  materials,  e.g.,  magnetostrictives,  piezoceramics,  and  Shape  Memory  Alloys  (SMAs),  hinders  the 
wider  applicability  of  such  materials  in  actuators.  A  fundamental  idea  in  coping  with  hysteresis  is  to  formulate  the 
mathematical  model  of  hysteresis  and  use  inverse  compensation  to  cancel  out  the  hysteretic  effect.  This  idea  can  be 
found  in  the  works  of  Tao  and  Kokotovic,1  Smith,2  and  Galinaitis  and  Rogers,3  to  name  a  few. 

Hysteresis  models  for  smart  materials  can  be  classified  into  physics  based  models  and  phenomenological  models. 
An  example  of  physics  based  model  is  the  Jiles- Atherton  model  for  ferromagnetic  hysteresis,4  where  hysteresis  is 
considered  to  arise  from  pinning  of  domain  walls  on  defect  sites.  The  most  popular  phenomenological  hysteresis 
model  used  in  control  of  smart  actuators  has  been  the  Preisach  model. 5-8  A  similar  type  of  operator,  called 
Krasnosel’skii-Pokrovskii  (KP)  operator  was  used  by  Galinaitis  and  Rogers  in  modeling  a  piezoelectric  actuator.3 
Though  the  Preisach  model  does  not  provide  any  physical  insight  into  the  problem,  it  provides  a  means  of  developing 
phenomenological  models  that  are  capable  of  producing  behaviors  similar  to  physical  systems  (see  Mayergoyz9  for 
an  excellent  exposition). 

In  this  paper  we  present  a  systematic  approach  for  control  of  smart  actuators.  The  method  is  illustrated  through 
the  example  of  controlling  a  commercially  available  magnetostrictive  actuator.  The  model  we  use  for  the  magne¬ 
tostrictive  actuator  is  based  on  earlier  work  of  our  group.10,11  It  was  shown  that  a  key  component  of  a  low-order 
model  is  as  shown  in  Figure  4.  The  rate-independent  hysteresis  operator  is  a  classical  Preisach  operator  followed  by 
a  squaring  operator.  The  input  of  the  hysteresis  operator  is  the  current  input  to  the  actuator,  and  the  output  of 
the  squaring  operator  is  a  quantity  with  the  dimension  of  force.  In  this  paper,  we  deal  with  low  frequency  inputs 
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of  less  than  5  Hz,  as  we  wish  to  focus  on  the  hysteresis  operator  only.  For  this  case,  the  linear  system  in  Figure  4 
approximates  to  a  constant  with  dimension  of  length/force. 

A  constrained  least-squares  algorithm  is  employed  to  identify  a  discrete  approximation  to  the  Preisach  measure 
(also  called  weighting  function  in  literature).  We  then  discuss  a  nonlinear  inversion  algorithm  for  the  resulting 
Preisach  operator,  based  on  the  theory  of  strictly-increasing  operators.  This  algorithm,  called  closest  match  algorithm, 
yields  a  control  input  signal  producing  the  closest  match  to  a  desired  magnetostrictive  response. 

Previous  work  closest  to  this  paper  was  by  Hughes  and  Wen.6-'  Although  both  their  papers  and  this  paper 
deal  with  the  themes  of  (a)  identification  of  a  Preisach  operator;  (b)  numerical  inverse  computation  of  the  Preisach 
operator  and  (c)  experimental  implementation,  there  are  significant  differences  both  in  the  method  and  complexity. 
First,  Hughes  and  Wen  use  first  order  reversal  curves  and  polynomial  fits  to  identify  a  continuous  Preisach  density 
function,  while  we  use  constrained  least  squares  method  to  identify  a  discrete  Preisach  measure.  Second,  Hughes 
and  Wen  use  an  implicit  function  theorem  approach  to  exactly  invert  the  Preisach  operator,  whereas  we  only  try  to 
find  the  input  yielding  the  closest  match  to  the  desired  output.  This  difference  is  significant  both  mathematically 
and  in  computational  savings.  In  Venkataraman  and  Krishnaprasad,12  we  have  shown  that  in  general,  one  can  only 
hope  to  approximately  invert  the  Preisach  operator  with  a  continuous  density  function,  and  trying  to  exactly  invert 
it  is  fraught  with  numerical  ill-conditioning.  The  implementation  of  the  closest  match  algorithm  is  both  numerically 
robust  and  time-saving. 

The  remainder  of  the  paper  is  organized  as  follows.  Section  2  provides  an  introduction  to  the  Preisach  operator, 
where  the  emphasis  is  on  how  to  evaluate  the  output  of  the  Preisach  operator  given  an  initial  memory  curve  and  the 
input.  Disretization  scheme  and  identification  algorithm  will  be  discussed  in  Section  3.  Section  4  is  devoted  to  the 
closest  match  inverse  algorithm,  which  fully  exploits  the  strictly-increasing  property  of  the  Preisach  operator  and  the 
discrete  structure  of  the  problem.  Experimental  results  are  given  in  Section  5.  Concluding  remarks  and  discussions 
on  possible  future  work  are  provided  in  Section  6. 

2.  THE  PREISACH  MODEL 

Consider  a  simple  hysteretic  element  (relay)  shown  in  Figure  1.  The  relationship  between  the  “input”  variable  u  and 
the  “output”  variable  v  at  each  instant  of  time  t  can  be  described  by: 

{v  =  +1  if  u>  a, 

v  =  —  1  if  u  <  /3,  .  (1) 

v  remains  unchanged  if  /3  <  u  <  a. 

Call  the  operator  relating  tt(-)  to  v(-)  as  7/3, «[«](•)>  where  we  now  view  the  input  and  output  variables  as  functions  of 
time.  This  operator  is  sometimes  referred  to  as  an  elementary  Preisach  hysteron  since  it  is  a  basic  block  from  which 
the  Preisach  operator  P[-]  will  be  constructed.  We  now  outline  this  construction.  Suppose  «(•)  €  C[0,  T]  is  the  input 
to  the  elementary  hysteron.  The  output  of  the  Preisach  operator  is  defined  as: 


u{t)  =T[u](t)  =  (  (  ii(fi,a)xfl p,a[u](t)dpda, 

J  J  ot>(3 


where  /x(-,  •)  is  a  density  function  (also  called  weighting  function  or  the  Preisach  measure).  This  representation  is  the 
most  natural  one  for  the  Preisach  operator  and  is  closest  to  Preisach’s  original  definition.9  The  Preisach  operator  has 
non-local  memory  and  it  “remembers”  the  dominant  maximum  and  minimum  values  of  the  past  input.  For  a  review 
of  this  and  other  basic  properties  of  the  Preisach  operator,  please  refer  to  Mayergoyz9  and  Brokate  and  Sprekels.13 

The  memory  effect  of  the  Preisach  operator  can  be  captured  by  curves  in  the  Preisach  (/3,  a)  plane.  To  simplify 
the  discussion,  we  assume  the  Preisach  measure  /z(-,  •)  has  a  compact  support,  i.e.,  a)  =  0  if  (3  <  fio  or  a  >  ao  for 

some  p0,  ao ■  Then  the  Preisach  plane  P  =  {(/3,  a)\a  >  ft,  ft  >  /?o,  a  <  ao},  as  shown  in  Figure  2(a).  Each  (/3,  a)  £  P 
is  identified  with  the  hysteron  7/3,a-  At  each  time  instant  t,  P  can  be  divided  into  two  regions: 


P_(t)  =  {(/?,  a)  £  P |  output  of  7^0,  at  t  is  —  1}, 
P+(t)  =  {(/?,  a)  £  P\  output  of  73jC,  at  t  is  +  1}, 


Figure  1.  Illustration  of  the  elementary  Preisach  hysteron 


so  that  P  =  PJt)  U  P+(t)  at  all  times.  Equation  (2)  can  be  rewritten  as: 

u(t)  =  f  f  a)d/3da  —  f  f  fi(/3,a)d/3da.  (3) 

J  Jp+(t)  J  Jp_(t) 

Now  assume  that  at  some  initial  time  to,  the  input  u(to)  =  Uq  <  /?o-  Then  the  output  of  every  hysteron  operator 
is  -1.  Therefore  PJt o)  =  P,  P+{to)  =  0  and  it  corresponds  to  the  “negative  saturation” (Figure  2(b)).  Next  we 
assume  that  the  input  is  monotonically  increased  to  some  maximum  value  at  t\  with  u(t\)  =  U\.  The  output  of 
is  switched  to  +1  as  the  input  u(t)  increases  past  a.  Thus  at  time  t\,  the  boundary  between  PJt i)  and  P+  (tj  is  the 
horizontal  line  a  =  rti(Figure  2(c)).  Next  we  assume  that  the  input  starts  to  decrease  monotonically  until  it  stops  at 
t2  with  u(t2)  =  u2.  It’s  easy  to  see  that  the  output  of  ~fp.a  becomes  -1  as  u(t)  sweeps  past  /?,  and  correspondingly, 
a  vertical  line  segment  (3  =  u2  is  generated  as  part  of  the  boundary  (Figure  2(d)).  Further  input  reversals  generate 
additional  horizontal  or  vertical  boundary  segments. 


a  a 


Figure  2.  Memory  curve  in  the  Preisach  plane 

From  the  above  illustration,  we  can  see  that  each  of  P_  and  P+  is  a  connected  set,  and  the  output  of  the  Preisach 
operator  is  determined  by  the  boundary  between  P_  and  P+-  The  boundary  is  also  called  the  memory  curve,  since  it 


provides  information  about  the  state  of  any  hysteron.  At  any  time  instant,  the  memory  curve  is  a  piecewise  constant, 
nonincreasing  function  of  (3.  We  also  note  that  due  to  its  staircase  structure,  the  memory  curve  is  fully  captured  by 
its  corner  points,  which  correspond  exactly  to  the  past  dominant  maximum  and  minimum  input  values.  Motivated 
by  this  observation,  we  store  and  update  only  these  corner  points  in  our  numerical  implementation  of  the  Preisach 
model. 


3.  IDENTIFICATION  OF  THE  PREISACH  MEASURE 
3.1.  Review  of  identification  methods 

A  classical  method  for  identifying  the  Preisach  measure  is  using  the  so  called  first  order  reversal  curves ,  detailed  in 
Mayergoyz.9  A  first  order  reversal  curve  can  be  generated  by  first  bringing  the  input  to  the  negative  saturation, 
followed  by  a  monotonic  increase  to  a,  then  a  monotonic  decrease  to  (3.  The  term  “first  order  reversal”  comes  from 
the  fact  that  each  of  these  curves  is  formed  after  the  first  reversal  of  the  input.  Denote  the  output  value  as  f([3,a) 
when  the  input  reaches  f3.  Then  the  measure  /x(/3,  a)  can  be  obtained  as 


fo  ^  1  d2f((3,a) 

^'a)  =  2-imr 


(4) 


Equatin  (4)  is  useful  only  when  the  two-dimensional  surface  f(/3,a)  is  twice  differentiable,  which  is  not  the  case 
for  measured  curves  in  experiments.  To  overcome  this  difficulty,  a  smooth  approximation  surface  f(/3,a)  is  fit  to 
the  data  points. 6-8  Hughes  and  Wen6-'  approximated  the  surface  by  polynomials  using  a  least  squares  method. 
Gorbet,  Wang  and  Morris  employed  functions  with  specific  forms,  and  the  parameters  were  obtained  via  a  weighted 
least  square  algorithm.8  As  pointed  out  in  Gorbet,  Wang  and  Morris,8  deriving  the  measure  by  differentiating  a 
fitted  surface  is  inherently  imprecise,  since  different  type  of  approximating  functions  lead  to  quite  different  measure 
distributions. 

Hoffmann  and  Sprekels14  proposed  a  scheme  to  identify  the  Preisach  measure  directly.  By  devising  the  input 
sequence  carefully,  they  set  up  independent  blocks  of  linear  equations  involving  the  output  measurements  and  the 
measure  masses  in  the  discretized  Preisach  plane.  Each  block  of  equations  can  be  solved  successively  to  obtain  the 
measure.  This  scheme  is  very  sensitive  to  experimental  errors  as  one  can  easily  see.  Using  the  identified  discrete 
measure,14  Hoffman  and  Meyer15  approximated  the  (continuous)  Preisach  measure  in  terms  of  a  set  of  basis  functions. 
A  least  squares  method  was  applied  to  compute  the  coefficients. 

Another  way  to  obtain  the  measure  is  driving  the  system  with  a  “reasonably”  rich  input  signal,  measuring  the 
output  and  then  estimating  the  Preisach  measure  by  a  least  squares  method.  This  idea  appeared  in  the  work  of 
Banks  and  his  colleagues,16'17  where  they  investigated  the  identification  problem  of  the  KP  operator.  Galinaitis 
and  Rogers3  used  the  same  idea  to  identify  the  weights  for  a  discretized  KP  operator.  We  will  also  adopt  the  least 
squares  method  for  measure  identificaiton  in  this  paper. 


3.2.  Identification  scheme 

Magnetostrictive  actuators,  due  to  the  capacity  of  the  windings  or  other  practical  reasons,  have  to  be  operated  with 
their  inputs  within  a  specific  range.  As  a  consequence,  we  will  not  be  able  to  visit  the  whole  Preisach  plane  and 
identify  the  measure  everywhere  during  the  identification  process.  We  assume  the  input  range  is  [umin,umax].  In 
Figure  3,  the  big  triangle  represents  the  whole  Preisach  plane,  while  the  smaller  triangle  is  the  region  we  can  visit 
and  we  denote  it  by  Hi.  The  region  outside  Hi  in  the  Preisach  plane  is  denoted  by  Ho-  Since  the  input  u(t)  never 
goes  beyond  the  limits,  states  of  hysterons  in  Ho  remain  unchanged.  Thus  the  bulk  contribution  to  the  output  from 
Ho  is  a  constant  and  we  denote  it  by  to o- 

The  input  is  discretized  into  L  +  1  levels  uniformly  and  we  label  the  cells  in  the  grid  as  illustrated  in  Figure  3 
for  L  =  9.  The  measure  mass  inside  each  cell  is  assumed  to  concentrate  at  the  cell  center.  The  quantities  we  want 
to  identify  include  measure  masses  iit]  ,i  =  1,  •  ■  • ,  L,  j  =  1,  •  •  • ,  i  and  loq.  To  simplify  the  discussion,  with  a  little  bit 
notation  abuse,  we  write  [Vij,i  =  1,  •  •  • ,  L,  j  =  1,  •  •  • ,  i  into  a  column  vector  fik,  k  =  1,  •  •  • ,  K,  where  K  =  rAL+'1'1 . 

To  initialize  the  states  of  hysterons,  we  first  increase  the  input  to  umax  and  then  reduce  it  to  umin.  This  sets 
the  state  of  each  hysteron  in  Hi  to  -1.  We  may  also  initialize  each  hysteron  in  Hi  to  +1  by  decreasing  the  input 
to  umin  followed  by  bringing  it  to  umax •  We  then  apply  some  piecewise  monotone  continuous  input  u(t),t  €  [0,T] 
which  contains  sufficiently  rich  information  (by  this  we  mean  Hi  should  be  visited  completely),  and  measure  the 


a 

t 


Figure  3.  Discretization  of  the  Preisach  plane  ( L  =  9) 

output  u)(t).  Signal  u(t),u>(t)  are  sampled  into  sequences  u[n],u)[n],n  =  1  ,  •  •  • , IV.  The  input  sequence  u[n]  (after 
discretization)  is  fed  into  the  discretized  Preisach  operator  and  the  state  of  each  hysteron,  7 k[n\,k  =  is 

computed.  The  output  of  the  Preisach  model  at  time  instant  n  can  be  expressed  as: 

K 

(5) 

the  parameters  are  determined  in  such  a  way 

(6) 

1 1  r  w[i]  ' 

1  u[2\ 

;  ,  b=  . 

1  J  [  w[IV]  _ 

Then  the  problem  becomes  finding  the  parameter  vector  9,  such  that 

II  A0-b\ | 

is  minimized,  where  ||  •  ||  stands  for  the  Euclidean  norm  in  1 ZN .  Since  we  require  /Xfc  >  0,  k  =  1,  •  •  • ,  K,  it  is  a  least 
square  error  optimization  problem  with  constraints. 

What  we  have  identified  is  a  discrete  approximation  to  the  Preisach  measure.  Sometimes  a  continuous  measure 
function  is  of  more  interest,  e.g.,  when  we  want  to  study  the  analytic  properties  of  the  system.  In  that  case,  we  may 
fit  a  smooth  function  from  the  identified  discrete  measure. 


urn 


= ^0 + y^M/c7fc[ 


k= 1 


where  [ik  s  are  yet  to  be  found. 

We  use  the  least  squares  method  to  estimate  the  parameters,  i.e. 
that 

N 


E 

n—  1 


\uj\n\  —  uj\n\ 


is  minimized.  This  can  be  written  in  a  more  compact  form.  Let 
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u>0 

7i[i]  72  [1] 
7i[2]  72  [2] 


7a- [1] 
7a:  [2] 


7i[iV]  72  [IV]  •••  7a- [IV] 


4.  INVERSION  OF  THE  PREISACH  OPERATOR 


The  general  structure  of  models  for  smart  actuators  that  capture  both  hysteresis  and  dynamic  behaviour  is  shown 
in  Figure  4  .  In  the  figure,  G(s)  represents  the  transfer  function  of  the  linear  part  in  the  actuator,  while  W  denotes 
a  rate-independent  hysteretic  nonlinearity.  Venkataraman10  has  shown  that  a  key  component  of  a  low-order  model 
for  magnetostriction  in  Terfenol-D  has  a  structure  resembling  Figure  4  . 


Rate-independent  Linear  system 

hysteresis  operator 


Figure  4.  Structure  of  models  for  smart  actuators 

A  basic  idea  for  controller  synthesis  for  such  systems  is  to  design  a  right  inverse  operator  W~1  for  W  as  shown 
in  Figure  5.  Then  w(-)  =  u(-)  and  the  controller  design  problem  is  reduced  to  designing  a  linear  controller  K(s)  for 
the  linear  system  G(s). 


y~ref  (■) 


Controller 


Figure  5.  Controller  design  schematic 

In  the  context  of  this  paper,  we  consider  W  as  a  Preisach  operator.  The  Preisach  operator  is  highly  nonlinear,  and 
in  general,  we  cannot  find  closed  form  formulas  for  the  inverse  operator,  unless  the  measure  function  is  of  some  special 
form,  like  in  the  work  of  Galinaitis  and  Rogers.18  Hughes  and  Wen6-'  utilized  the  first  order  reversal  curves  /(/? ,  a) 
in  computing  the  numerical  inverse  of  the  Preisach  operator.  They  defined  F(/3,  a)  =  |[/(a,  a)  —  f(/3,  a)]  (called  the 
Everett  surface).  For  output  change  A w, 

Aw  =  2  F(/3,a)  (7) 

defines  implicitly  two  inverse  maps  Ga(-,  •)  and  G0(-,  •): 

at.  =  Ga(Aw,  (3),  (8) 

f3  =  G0(/3,Ato).  (9) 

Ga(-,  •)  and  Gp(-,  •)  is  then  used  to  compute  the  desired  input  given  the  desired  output  and  past  input  history.  This 
method  relies  on  measurement  of  all  first  order  reversal  curves  and  therefore  is  subject  to  experimental  errors.  Also 
computing  Ga ,  G0  involves  solving  nonlinear  equations  and  therefore  is  not  a  trivial  task. 

4.1.  Contraction  mapping  algorithm  for  inversion 

Under  mild  assumptions,  the  Preisach  operator  is  Lipshitz  continuous  and  incrementally  strictly  increasing.13  Venkatara¬ 
man  and  Krishnaprasad1912  exploited  these  properties  and  proposed  an  inversion  algorithm  for  the  Preisach  operator 
based  on  the  contraction  mapping  principle.  The  algorithm  is  summarized  in  the  following  theorem: 


Theorem  4.1: [Contraction  Mapping  Algorithm  for  Inversion]1912  Let  X  =  Ci[0,T\,  wherel  =  [a,b\. 
Ci[0,T]  denotes  the  space  of  continous  functions  defined  on  [0, T]  taking  values  in  I.  Let  F  :  X  — »•  Y,  be  an 
incrementally  strictly  increasing,  strongly  Lipshitz  continuous  Preisach  operator  (with  Lipshitz  constant  k2)  with 
some  initial  memory  curve  ip~\,  and  Y  is  the  range  ofT.  Given  e  >  0  and  the  operator  equation 

T(x)  =  y,  (10) 


where  y  £  Y  is  known,  consider  the  algorithm: 

•  pick  any  Xo  £  X; 

•  while  ||  xn  —  xn-\  ||  >  e: 


xn+i  =  z-(y  +  (k2xn  -  r(xn))). 
k2 


(11) 


The  sequence  {xn}  terminates  at  z  which  satisfies  |j  z  —  x  ||<  e  where  x  is  the  solution  of  (10).  The  rate  of  convergence 
is  linear. 

When  the  Preisach  measure  is  discretized,  the  Preisach  operator  is  no  longer  Lipschitz  continuous  and  therefore 
the  contraction  mapping  algorithm  for  inversion  does  not  work  efficiently.  Indeed,  since  the  output  of  P  at  any  time 
instant  can  take  only  finite  number  of  possible  values  (due  to  the  finite  number  of  measure  atoms),  from  (11),  we  see 
that  for  almost  all  y, 

II  xn  -  xn—\  ||=  —  II  y-  r(xn_i)  II >Sy>  0. 

In  addition,  the  discrete  measure  masses  are,  in  general,  not  uniform.  These  factors  make  it  difficult  to  choose  an 
appropriate  stopping  criterion  e:  picking  e  big  we  lose  accuracy;  picking  e  small  we  will  get  stuck  if  e  <  Sy.  Note 
that  these  remarks  don’t  invalidate  Theorem  4.1:  Theorem  4.1  works  perfectly  for  continuous  measure  distribution. 
However,  we  need  a  more  practical  inversion  algorithm  for  the  discrete  measure  case. 


4.2.  Closest  match  algorithm  for  inversion 

We  propose  a  new  inversion  algorithm  in  this  subsection.  This  algorithm,  like  the  contraction  mapping  algorithm,  is 
also  based  on  the  strictly  increasing  property  of  the  Preisach  operator.  It  fully  utilizes  the  discrete  structure  of  the 
problem.  We  name  it  closest  match  algorithm  because  it  always  generates  input  whose  output  matches  the  desired 
output  most  closely  among  all  possible  inputs. 

Note  that  due  to  discretization,  the  input  can  only  take  values  from  a  finite  set  U  =  {ui,  1  <  l  <  L  +  1}  with 
each  ui,  1  <  l  <  L  +  1,  representing  an  input  level.  To  be  precise,  let 

a  nmax  nmin 


then  ui  =  umin  +  {l  —  1)  A  u.  Thus  tti  =  umin  and  ul+ i  =  umax-  The  inversion  problem  is:  given  an  initial  memory 
curve  ^>(°)(ffom  which  the  initial  input  and  output  u/°)  can  be  derived)  and  a  desired  output  uj,  find  u*  £  U, 
such  that 

|P(u*;  ip^)  —  w\=  min  |F(u;  ip^)  —  u>\.  (12) 

u£U 

Also  the  algorithm  should  return  the  resulting  memory  curve  ip*  for  later  use.  Note  in  (12)  we  explictly  put  ip as 
argument  of  P  to  emphasize  the  effect  of  the  memory  curve  on  the  output. 

The  intuitive  idea  is  as  follows.  Consider  the  case  where  the  current  output  is  less  than  the  desired  value  u> 
(the  case  >  u>  is  treated  in  exactly  the  same  way  with  some  obvious  modification  ).  We  keep  increasing  the  input 
by  one  level  in  each  iteration  until,  say  at  iteration  n,  the  input  uljli  reaches  umax ,  or  the  output  corresponding 
to  u exceeds  uj.  For  the  first  case,  the  optimal  input  is  clearly  umax',  for  the  second  case,  two  candidates  for  the 
optimal  input  u*  are  u^n~^  and  .  We  then  take  u*  to  be  the  one  with  smaller  output  error.  Note  that  we  need 
back  up  the  memory  curve  whenever  we  increase  input,  so  that  we  can  always  retrieve  the  consistent  memory  curve 
with  u* . 

We  now  describe  the  algorithm  in  detail. 


Closest  Match  Algorithm: 


•  Step  0  [Initialization].  Set  n  =  0.  Compare  u/°)  and  co:  if  co ^  =  co,  let  u*  =  u go  to  Step  3;  if 

<  Co,  go  to  Step  1.1;  otherwise  go  to  Step  2.1; 

•  Step  1  [Case  <  Co]. 

-  Step  1.1:  If  u ^  =  umax,  let  u*  =  u^n\ip*  =  ip(n\  go  to  Step  3;  otherwise  u^n+1^  =  +  A tt,V>  = 

^l"l[back  up  the  memory  curve],  n  =  n  +  1,  go  to  Step  1.2; 

—  Step  1.2:  Evaluate  =  r(it!ra); ^Mn_1)),  and  (at  the  same  time)  update  the  memory  curve  to  if>(n\ 
Compare  with  a):  if  a/")  =  Co ,  let  u*  =  u^n\ip*  =  ip( n\  go  to  Step  3;  if  co^  <  Co,  go  to  Step  1.1; 
otherwise  go  to  Step  1.3; 

-  Step  1.3:  If  | —  cZ> |  <  —  Co\,  let  u*  =  u^n\tj)*  =  ip(n'1 ,  go  to  Step  3;  otherwise  u*  = 

[restore  the  memory  curve],  go  to  Step  3; 

•  Step  2  [Case  >  Co]. 

—  Step  2.1:  If  =  umin,  let  u*  =  u^n\%p*  =  ip^n\  go  to  Step  3;  otherwise  —  A u, ?/>  = 

</)l"l[back  up  the  memory  curve],  n  =  n  +  1,  go  to  Step  2.2; 

—  Step  2.2:  Evaluate  a/ral  =  r(tt!n); ^Mn_1l),  and  (at  the  same  time)  update  the  memory  curve  to  %p(n\ 
Compare  lo^  with  Co:  if  =  Co,  let  u*  =  u(n\i/j*  =  go  to  Step  3;  if  to^71'1  >  Co,  go  to  Step  2.1; 
otherwise  go  to  Step  2.3; 

—  Step  2.3:  If  |u/ra)  —  Co\  <  |u/n-1)  —Co\,  let  u*  =  u (n\ip*  =  ip(n') ,  go  to  Step  3;  otherwise  u*  =  u(n~1'>,ip*  = 
[restore  the  memory  curve],  go  to  Step  3; 

•  Step  3.  Exit. 

It’s  not  hard  to  see  the  above  algorithm  yields  the  best  input  u*  in  at  most  L  iterations.  And  in  each  iteration, 
the  evaluation  of  to ^  is  very  fast  since  the  input  has  changed  by  one  level  and  thus  we  need  only  update  states  of 
hysterons  corresponding  to  that  level.  In  other  words,  the  state  of  each  hysteron  needs  to  be  updated  at  most  once 
during  the  whole  process  of  finding  u*.  These  factors  combine  to  make  this  algorithm  simple  and  efficient. 

Other  types  of  search  algorithm,  like  bisection  algorithm  may  also  be  applied  to  find  u*.  However,  bisection 
algorithm  may  involve  re-evaluation  of  states  of  some  hysterons  for  many  times,  which  makes  the  algorithm  slow. 

The  closest  match  algorithm  generates  input  in  a  discrete  value  set  and  thus  apparently  it  is  discontinuous  in  time. 
This  is  not  true  because  we  tacitly  assume  the  input  is  changed  continuously  and  monotonically  in  each  iteration  in 
evaluating  the  output  of  the  discretized  Preisach  operator.  This  continuity  is  achieved  in  real  system  implementation 
by  linearly  interpolating  between  the  computed  values. 

5.  EXPERIMENTAL  RESULTS 

In  this  section,  we  will  apply  the  identification  and  inversion  schemes  to  open  loop  control  of  a  magnetostrictive 
actuator.  Magnetostriction  is  the  phenomenon  of  strong  coupling  between  magnetic  properties  and  mechanical 
properties  of  some  ferromagnetic  materials  (e.g.,  Terfenol-D):  strains  are  generated  in  response  to  an  applied  magnetic 
field,  while  conversely,  mechanical  stresses  in  the  materials  produce  measurable  changes  in  magnetization.  Figure  6 
shows  a  sectional  view  of  a  Terfenol-D  actuator  manufactured  by  ETREMA  Products,  Inc.  By  varying  the  current 
in  the  coil,  we  vary  the  magnetic  field  in  the  Terfenol-D  rod  and  thus  control  the  motion  of  the  rod  head.  Figure  7 
displays  the  hysteresis  in  the  magnetostrictive  actuator. 

As  mentioned  in  Section  1,  the  magnetostriction  A  is  connected  to  the  magnetization  M  by 

A  =  a\M2.  (13) 

We  can  identify  the  coefficient  a\  from  the  saturation  magnetization  Ms  and  the  saturation  magnetostriction  As 
provided  by  the  manufacturer: 

As 
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Figure  6.  Sectional  view  of  the  Terfenol-D  actuator 


Figure  7.  Hysteresis  in  the  magnetostrictive  actuator 


Let  the  length  of  the  magnetostrictive  rod  be  lrod-  Given  a  measurement  of  the  displacement  d,  since  d  =  lrod A, 
the  underlying  magnetization  M  is  determined  by, 

M  =  ±\f^7d'  (14) 

and  the  sign  of  M  is  determined  with  further  information  on  the  input.  The  applied  magnetic  field  H  is  related  to 
the  input  current  through 

H  =  NI  +  Hbias,  (15) 

where  N  is  the  number  of  coils  per  unit  length  and  H^ias  is  the  bias  magnetic  field  produced  by  permanent  magnets. 
Huas  is  necessary  for  generating  bidirectional  strains  and  it  can  be  identified  easily. 

We  will  treat  the  magnetic  field  H  as  input  and  the  bulk  magnetization  M  as  output  using  transformations  (14) 
and  (15).  And  we  employ  the  Preisach  operator  to  model  the  hysteretic  relationship  between  M  and  H.  We  will 
identify  the  Preisach  measure,  and  then  carry  out  inverse  compensation  based  on  the  identified  measure.  An  open 
loop  tracking  experiment  will  be  done  to  check  the  performance  of  the  identification  and  inversion  algorithms. 

5.1.  Experiment  setup 

Our  experimental  setup  is  as  shown  in  Figure  8  .  DSpace  ControlDesk  is  a  powerful  tool  for  real-time  simulation  and 
control.  It  can  generate  system  models  from  Simulink  of  Matlab,  download  real-time  applications  into  a  DSP  board, 
monitor  and  control  a  system  in  real  time  by  collecting  data  and  sending  out  commands.  The  data  it  collects  can 
be  displayed  or  be  saved  on  disk  for  post-processing.  The  displacement  of  the  actuator  is  measured  with  a  LVDT 


sensor. 
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Figure  8.  Experiment  setup 


Figure  9.  Distribution  of  the  identified  measure 

5.2.  Measure  identification  and  validation 

The  magnetic  field  input  H  is  limited  in  the  range  [-40  Oe,  480  Oe]  and  is  discretized  into  26  levels.  Figure  9  shows 
distribution  of  the  identified  measure.  The  constant  contribution  from  Do  (see  Figure  3)  is  estimated  to  be  0.3466. 

To  verify  the  identified  measure,  we  apply  same  input  signal  (which  differs  from  the  input  used  for  identification) 
to  the  actuator  and  the  Preisach  model.  Figure  10  shows  the  comparison  of  the  actuator  output  and  the  output  of 
the  Preisach  model.  We  can  see  they  agree  reasonably  well  and  therefore  the  identified  Preisach  operator  provides  a 
good  approximation  to  the  actuator. 

5.3.  Open  loop  tracking  experiment 

Finally,  we  do  an  open  loop  tracking  experiment  to  test  the  overall  performance  of  our  identificaiton  and  inversion 
scheme.  Given  a  desired  trajectory,  whose  amplitude  and  frequency  are  both  varying,  we  compute  a  desired  input 
using  the  closest  match  inversion  algorithm.  The  computed  input  signal  is  then  sent  to  the  actuator  and  the 
displacement  trajectory  is  measured.  Figure  11  shows  the  comparison  of  the  desired  trajectory  and  the  actual 
trajectory.  As  we  can  see,  the  tracking  error  lies  within  a  small  interval  [— 3/r?n,  3/rm]. 


6.  CONCLUSIONS  AND  FUTURE  WORK 

In  this  paper,  we  have  proposed  a  systematic  approach  for  control  of  smart  actuators.  We  model  the  hysteresis  by 
the  Preisach  operator  and  identify  the  Preisach  measure  using  a  constrained  least  square  method.  We  have  presented 
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Figure  10.  Validation  of  the  identified  measure 
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Figure  11.  Open  loop  tracking  result 


an  inversion  algorithm  based  on  the  strictly  increasing  property  of  the  Preisach  operator.  The  effectiveness  of  our 

approach  is  demonstrated  via  an  open-loop  trajectory  tracking  experiment. 

There  are  a  couple  of  possible  and  interesting  directions  to  extend  this  work: 

•  It  is  well  known  that  properties  of  smart  actuators  may  vary  with  time,  temperature,  etc.  This  means  we  might 
need  to  re-identify  the  model  quite  often  or  even  on-line  identification  is  necessary.  Also  Off-the-shelf  algorithms 
for  solving  the  constrained  least  squares  problem  can  be  very  time-consuming  when  the  discretization  gets  fine. 
Therefore  a  fast  and  efficient  identification  algorithm  will  be  very  useful. 

•  In  this  paper  we  consider  low  frequency  input  signals  and  thus  ignore  the  dynamic  part  G(s)  in  the  model 
(Figure  4).  Extending  of  current  work  to  accomodate  wider  frequency  bandwidth  is  of  practical  importance. 

REFERENCES 

1.  G.  Tao  and  P.V.  Kokotovic,  “Adaptive  control  of  plants  with  unkown  hystereses,”  IEEE.  Trans.  Automat.  Contr. 
40(2),  pp.  200-212,  1995. 

2.  R.C.  Smith,  “Inverse  compensation  for  hysteresis  in  magnetostrictive  transducers,”  CRSC  Technical  Report, 
North  Carolina  State  University  ,  CRSC-TR98-36,  1998. 

3.  W.S.  Galinaitis  and  R.C.  Rogers,  “Control  of  a  hysteretic  actuator  using  inverse  hysteresis  compensation,”  in 
Mathematics  and  Control  in  Smart  Structures,  V.  Varadan,  ed.,  Proc.  SPIE  3323,  pp.  267-277,  1998. 

4.  D.C.  Jiles  and  D.L.  Atherton,  “Theory  of  ferromagnetic  hysteresis,”  J.  Magn.  Magn.  Mater.  61,  pp.  48-60, 
1986. 

5.  A. A.  Adly,  I.D.  Mayeygoyz,  and  A.  Bergqvist,  “Preisach  modeling  of  magnetostrictive  hysteresis,”  J.  Appl. 
Phys.  69(8),  pp.  5777-5779,  1991. 

6.  D.  Hughes  and  J.T.  Wen,  “Preisach  modeling  and  compensation  for  smart  material  hysteresis,”  in  Active  Ma¬ 
terials  and  Smart  Structures,  Proc.  SPIE  2427,  pp.  50-64,  1994. 

7.  D.  Hughes  and  J.T.  Wen,  “Preisach  modeling  of  piezoceramic  hysteresis;  independent  stress  effect,”  in  Mathe¬ 
matics  and  Control  in  Smart  Structures,  V.  Varadan,  ed.,  Proc.  SPIE  2442,  pp.  328-336,  1995. 

8.  R.B.  Gorbet,  D.W.L.  Wang,  and  K.A.  Morris,  “Preisach  model  identification  of  a  two-wire  SMA  actuator,”  in 
Proceedings  of  IEEE  International  Conference  on  Robotics  and  Automation,  pp.  2161-2167,  1998. 

9.  I.D.  Mayergoyz,  Mathematical  Models  of  Hysteresis,  Springer  Verlag,  New  York,  1991. 

10.  R.  Venkataraman,  Modeling  and  Adaptive  Control  of  Magnetostrictive  Actuators,  Ph.  D.  thesis,  University 
of  Maryland,  College  Park,  available  at:  http  :  // www.isr .umd.edu/Tech.Reports / 1 S R/1999 / PhD  99  — 
1/ PhD  99  -  l.phtml,  1999. 

11.  R.  Venkataraman  and  P.S.  Krishnaprasad,  “A  model  for  a  thin  magnetostrictive  actuator,”  in  Proceedings  of 
the  Conference  on  Information  Sciences  and  Systems,  1998.  Also  published  as  a  technical  report,  TR  98-37, 
Institute  for  Systems  Research,  University  of  Maryland  at  College  Park. 

12.  R.  Venkataraman  and  P.S.  Krishnaprasad,  “Approximate  inversion  of  hysteresis:  theory  and  numerical  results,” 
in  Proceedings  of  the  39th  IEEE  Conference  on  Decision  and  Control,  pp.  4448-4454,  2000. 

13.  M.  Brokate  and  J.  Sprekels,  Hysteresis  and  Phase  Transitions,  Springer  Verlag,  New  York,  1996. 

14.  K.-H.  Hoffmann  and  J.  Sprekels,  “Identification  of  hysteresis  loops,”  J.  Comput.  Phys.  78,  pp.  215-230,  1988. 

15.  K.-H.  Hoffmann  and  G.H.  Meyer,  “A  least  squares  method  for  finding  the  Preisach  hysteresis  operator  from 
measurements,”  Numer.  Math.  55,  pp.  695-710,  1989. 

16.  H.T.  Banks,  A.J.  Kurdila,  and  G.  Webb,  “Identification  of  hysteretic  control  influence  operators  representing 
smart  actuators,  part  I:  formulation,”  Mathematical  Problems  in  Engineering  3(4),  pp.  287-328,  1997. 

17.  H.T.  Banks,  A.J.  Kurdila,  and  G.  Webb,  “Identification  of  hysteretic  control  influence  operators  representing 
smart  actuators,  part  II:  convergent  approximations,”  Journal  of  Intelligent  Material  Systems  and  Structures 
8(6),  pp.  536-550,  1997. 

18.  W.S.  Galinaitis  and  R.C.  Rogers,  “Compensation  for  hysteresis  using  bivariate  Preisach  models,”  in  Mathematics 
and  Control  in  Smart  Structures,  V.  Varadan  and  J.  Chandra,  eds.,  Proc.  SPIE  3039,  pp.  538-547,  1997. 

19.  R.  Venkataraman  and  P.S.  Krishnaprasad,  “A  novel  algorithm  for  the  inversion  of  the  Preisach  operator,”  in 
Mathematics  and  Control  in  Smart  Structures,  V.  Varadan,  ed.,  Proc.  SPIE  3984,  pp.  404-414,  2000. 


